Mobility of high-power solitons in saturable nonlinear photonic lattices 
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We theoretically study the properties of one-dimensional nonlinear saturable photonic lattices exhibiting 
multiple mobility windows for stationary solutions. The effective energy barrier decreases to a minimum in 
those power regions where a new intermediate stationary solution appears. As an application, we investigate 
the dynamics of high-power gaussian-like beams finding several regions where the light transport is enhanced. 
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For more than a decade, nonlinear waveguide arrays 
(WAs) have become an excellent experimental scenario, 
where to study the phenomenology appearing in periodic 
and non-periodic nonlinear dynamical systems [1-3]. Dif- 
ferent geometries, dimensions, and nonlinearities have 
been studied showing interesting and very different phe- 
nomenologies compared to continuous systems. For WAs 
with Kerr nonlinearity, a coupled mode approach leads 
to a discrete nonlinear Schrodinger (DNLS) equation. 
For photorefractive media with saturable nonlinearity 
the corresponding model is a s-DNLS [4]. 

Mobility of localized solutions in nonlinear cubic WAs 
is well known. As far as the power stays low, the energy 
barrier imposed by the discreteness and the nonlinearity 
(usually called Peierls-Nabarro (PN) potential [5]) will 
stay small and solutions will move across the lattice by 
just giving them a judicious kick [6,7]. For larger powers, 
the energy barrier grows and mobility is not possible any- 
more. However, for saturable systems there are several 
regions of power where the energy difference between the 
two fundamental localized solutions - the one centered at 
one site (odd mode) and the one centered between two 
sites (even mode) - vanishes for different power values [4] . 
Close to these points, there are regions of stability ex- 
change between the even and odd solutions. Since those 
regions exhibit bistability, the appearance of an inter- 
mediate asymmetric and unstable solution is inevitable, 
as it was shown in Ref. [8, 9] for two-dimensional sys- 
tems. Therefore, in this case, the effective energy barrier 
will strongly depend on the intermediate solutions (IS). 
For saturable one-dimensional (ID) systems, this cru- 
cial issue has not been clearly identified yet. We believe 
this element could be one of the keys to experimentally 
observe, to the best of our knowledge for the first time, 
good soliton mobility in ID nonlinear saturable photonic 
lattices. 

It was shown in Ref. [4] that the PN barrier be- 
comes minimal exactly at the points where the funda- 



mental solution's energies coincide. Contrary to what is 
expected, the fundamental solutions remain immobile in 
these points. We demonstrate that in order to achieve 
a good mobility it is necessary to increase the amount 
of power up to the bifurcation point where the IS dis- 
appears. A constraint method [9, 10] is used to identify 
the ISs and describe a pseudo-potential landscape among 
all stationary modes. By using the system properties, we 
found recurrent resonant behavior in power for gaussian- 
like shaped pulses showing enhanced mobility. 
In a ID system the s-DNLS is given by 

l-T~ + (Un+l + Un-1) ~ 12 = ' W 

dz 1 + \u n \ z 

where u n represents the light amplitude at site n, 7 the 
strength of the nonlinearity with respect to the cou- 
pling coefficient, and z a normalized propagation dis- 
tance along the waveguides. In order to understand the 
main phenomenology of these lattices, we first look for 
stationary solutions of the form u n (z) = u n exp(i\z), 
where u n G R and A is the propagation constant 
or frequency. Small-amplitude plane waves define the 
band A G [—2 — 7,2 — 7], while high-amplitude plane 
waves define a second band A G [—2,2]. Therefore in- 
phase stationary localized solutions are limited to ex- 
ist in the region A G [2 — 7,2], bifurcating from the 
fundamental modes of those bands [8]. Model (1) has 
two dynamically conserved quantities, the Hamiltonian 

H = ~ [En(Vi< +w n <+i) -7log(l + Kl 2 )] and 
the optical power P = \u n \ 2 . 

In the following, we discuss these properties for 7 = 10 
(focusing nonlinearity). A different 7 will produce differ- 
ent curves but the main saturable phenomenology will be 
preserved. Localized solutions are computed by using a 
standard Newton- Raphson method. Fig.l shows a power 
versus frequency diagram for both fundamental modes - 
the odd and the even solutions - including the IS. As 
expected, the IS corresponds to a non-symmetric profile 
connecting the two fundamental modes [see Fig.l-inset]. 
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Fig. 1. (Color online) P versus A for odd, even, and in- 
termediate solutions in full, dashed, and thin lines, re- 
spectively. Inset: profiles corresponding to filled circles. 



The two fundamental solutions cross each other - repeat- 
edly - as P increases. In regions where \P dd — Peven\ 
is large, a family of IS appears. It is remarkable that 
all the ISs of the first family have A = and connect 
both (even and odd) modes for this value. In such a sit- 
uation, the stationary solutions of Eq.(l) coincide with 
the ones of the integrable Ablowitz-Ladik equation [11], 
which has an analytic mobile solution. However, in the 
physical s-DNLS model, it is not expected to find radia- 
tionless travelling solutions, since mobile modes need to 
have the same power and not the same frequency. We 
perform a standard linear stability analysis [12] by com- 
puting the largest imaginary part ( "g" ) of the eigenvalue 
spectrum. Fig. 2(a) shows our results where g = implies 
stable solutions and g > unstable ones. For all regions 
where the fundamental solutions are simultaneously sta- 
ble (three regions in this plot), the unstable IS appears. 
In these regions there are points where the energy of both 
fundamental solutions is exactly the same [see inset in 
Fig. 2(a) for the first bistable region (P ~ 20)]. However, 
the effective energy barrier is not zero if the IS is consid- 
ered as well. In Fig. 2(b) we plot AHo = \H Q dd — Heven\ 
and AH = \H max — H min \ versus power. For the first two 
"bistable" regions, we clearly see that AHq goes to zero 
as it was previously predicted in Ref. [4]. However, there 
is always a nonzero barrier (AH) for the solution, which 
can be very small but it is - strictly speaking - nonzero. 
A first guess could be, that the most favorable region for 
mobility would be the one where AH is a minima. How- 
ever, this is not the case for stationary solutions. If we 
kick an odd or even mode, we are putting in motion an 
immobile-defined solution, therefore there is always radi- 
ation from tails. As a consequence, the power of the mov- 
ing solution is lower than the initial one. So, if we initially 
take the solution where AH is a minima (P m ), the effec- 
tive barrier will increase [see Fig. 2(b)]. A better option 
would be the one where P > P m where, due to radiation, 
the power and the effective barrier decrease. Now, in or- 
der to go deeper in the understanding of the dynamics 
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Fig. 2. (Color online) (a) g versus P for odd, even, and 
intermediate solutions in full, dashed, and thin lines, re- 
spectively. Inset: H versus P. (b) AHq (thick line) and 
AH (thin line) versus P. (c) Energy surface where light 
(dark) color denotes a high (low) i^-value. 



of ID saturable WAs, we construct an energy landscape. 
By defining the center of mass X = ^2 n n\u n \ 2 / P and 
using a constraint method [9,10], we compute H versus 
X and P. Fig. 2(c) shows our computations around the 
first bistable region. In this plot, X = n and n + 1 cor- 
respond to odd modes while X = n + 0.5 corresponds 
to even one. For low power, the potential is cubic-like in 
the sense that the odd mode is stable while the even one 
is unstable (lower curve for P = 15). By increasing the 
power, both fundamental solutions are simultaneously 
stable because they are both a local minima in this po- 
tential. Consequently, a maxima in-between appears, the 
unstable IS (middle curve for P = 20.5 where AHq « 0). 
Then, by further increasing the power, the odd solution 
transforms into an unstable maxima while the even one 
becomes the only minima of this potential (upper curve 
for P = 26). The IS originates when the even mode sta- 
bilizes; then it changes its center of mass to an odd mode 
(symmetrically to the right and to the left due to the sys- 
tem symmetry) . Finally, the IS disappears when the odd 
mode destabilizes. 

To the best of our knowledge, mobility in this kind 
of systems was never predicted for a more realistic ex- 
perimental input condition like gaussian input profiles. 
Previous simulations, starting from stationary solutions, 
observed good mobility [4,8,9,11]. However, saturable 
solutions are not well localized in the power-exchange 
regions and, furthermore, by increasing the power they 
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become broader. Therefore, in a experiment, dynamics 
will be strongly determined by the power and shape of 
the chosen beam profile. We took as an input beam a five- 
site wide gaussian-like profile: u n (0) = Aexp[— a(n — 
n c ) 2 ] exp[ik(n — n c )] for n = n Cl n c db l,n c db 2 and 
^n(O) = 0, otherwise. The kick k is proportional to the 
experimental angle and it does not alter the power but 
adds a small amount of effective kinetic energy. With 
this initial condition, we numerically integrate model (1) 
from z = 0toz = Zf and measure the center of mass 
at the lattice output: Xf = X(zf). Fig. 3(a) shows our 
results for different input power (P). For very low power 
(P ~ 0), the system is linear and mobility decreases 
as system becomes nonlinear up to P ~ 50, when the 
saturable system behaves as a cubic one [see Fig. 3(b)]. 
Then, by further increasing the power, solutions start 
to move. When the power is increased, fundamental so- 
lutions are geometrically similar and differences in the 
Hamiltonian are very small, therefore the profile is al- 
lowed to move with k ^ 0, as it is observed in the av- 
erage tendency of the curve [diagonal straight line in 
Fig. 3(a)]. However, some "resonant" dynamics is found 
for different levels of power. There are different regions 
where mobility is enhanced as shown in Figs. 3(c) and 
(d), with almost undamped motion. It is plausible to as- 
sume that this behavior corresponds to a manifestation 
of the continuously repeating bistable regions discussed 
for stationary solutions. Therefore, for even higher pow- 
ers, there should be good mobility, in principle without 
limitations. The mobility windows for gaussian pulses 
do not coincide with the mobility regions for station- 
ary solutions, since a gaussian pulse will always have a 
higher cost in power loss due to radiation and due to its 
shape (gaussian profiles do not match in shape and power 
with any fixed stationary solution, what implies different 
level of powers to observe similar dynamics). But, never- 
theless, the recurrent appearance of enhanced mobility 
for certain powers shows an excellent phenomenological 
agreement between gaussian and stationary profiles. 

To conclude, we have shown that in nonlinear sat- 
urable ID photonic lattices there are several regions of 
bistability where stationary solutions possess a small 
but nonzero energy barrier. The effective energy barrier 
among all stationary localized solutions was constructed 
allowing us to get a deeper understanding of discrete sat- 
urable nonlinear systems. By using these properties with 
a more realistic input condition, we were able to observe 
very good mobility and also to find different regions of 
resonant response where the mobility is enhanced. We 
hope that these findings will motivate experimentalists 
to explore this direction of high-power mobility, which 
is currently a drawback for the implementation of non- 
linear WAs in realistic all-optical operations in different 
power regimes. 
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Fig. 3. (Color online) (a) Output Xf versus input P. 
(b)-(d) Dynamical examples for P = 54, 300, and 550, 
respectively. 7 = 10, n c = 26, k = 0.3, a = 1/3, Zf = 50. 
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